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By analyzing simple models of fermions in lattice potentials we argue that the zero-temperature 
pairing instability of any ideal band-insulator occurs at a finite momentum. The resulting supersolid 
state is known as "pair density wave". The pairing momentum at the onset of instability is generally 
incommensurate as a result of phase-space restrictions and relative strengths of interband and in- 
traband pairing. However, commensurate pairing occurs in the strong-coupling limit and becomes a 
Cooper-channel analogue of the Halperin-Rice exciton condensation instability in indirect bandgap 
semiconductors. The exceptional sensitivity of incommensurate pairing to quantum fluctuations can 
lead to a strongly-correlated insulating regime and a non-BCS transition, even in the case of weak 
coupling as shown by an exact renormalization group analysis. 

I. INTRODUCTION 

One of the central challenges in the theory of high temperature superconductors is how to reconcile the mean- 
I I field character of Cooper pairing which is pivotal in any fermionic superfluid, and the seemingly fluctuation-driven 



f^ • superconducting transition in underdoped cuprates- . Motivated by this problem, we consider s-wave pairing instability 
^ in generic models of band-insulators and argue that certain phenomena familiar from cuprates, such as density-wave 
C/3 . patterns^ and fluctuation-driven transitions, are possible in these simple systems. The most direct physical realization 

of our models can be found in ultra-cold gases of alkali atoms tuned near a broad Feshbach resonance^. 

The simplest generic band-insulator is obtained at zero temperature by placing an even number of fermionic particles 

in each unit-cell of a periodic potential Vr '■ 



ci ^ -^ - /^ + ^r j c, - Uc\clc^c^ 



(1) 



^ ! H = I d^r 

We require that the chemical potential fj, sit in a band-gap, but allow changes of fi, lattice amplitude V and attractive 
interaction strength U to drive phase transitions (we use h = 1). In order to stay focused on universal aspects of 
^ ' unconventional pairing, we choose to work in the unitarity limit which defines a strongly-interacting quantum critical 
^T^ ■ point at which the only relevant interaction, included in ([T]), is related to the two-body vacuum scattering length a 
through a cut-off scale implicit in the following integral: 



1 m Y^ f d k 1 

U ^ "i^ ^ ^ J (27r)3 2e„fc v=o ' ^' 

n 

Here n and k are band-index and crystal momentum quantum numbers of quasiparticles, and tnk are their bare 
energies. 



X 

\-^ ■ II. PAIR DENSITY WAVE INSTABILITY 

In our recent work^ we argued that the pairing instability in the continuum model ([l} generally occurs at a finite 
incommensurate wavevector Q. Before subjecting this conclusion to more scrutiny, we outline the robust reasons for 
finite-momentum pairing. The mean-field inverse pairing susceptibility matrix 11 at frequency f2 = in the insulating 
state is: 



T11T12 



where G are reciprocal lattice vectors, ^nk = £nk — fJ-, and /(^) is Fermi-Dirac distribution function. The conservation 
of crystal momentum implies ^cq.G'q' = ^gg'{q) ^ (27i')^(5(<7 — q'). This expression is obtained by integrating out 
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FIG. 1: (a) n(q) in d = 3 for the potential Vr = V"^cos{2'jTXi/aL); Er='K^h^ /2ma\, Nb is the number of bands, No is the 

1=1 
number of order parameters, (b) A d = 2 phase diagram at T = 0. 

the internal Matsubara frequency in the bubble Feynman diagram, where the fermion propagators are written in the 
basis of Bloch states^. Such a basis choice introduces non-trivial momentum and band dependent vertex functions: 



niki:n2k2 
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-i(q+G)7- 



V'nifei ir)i^n2k2 i^) « ^C^l + ^2 - q) 



(4) 



Typical plots of the lowest eigenvalue n(q) of ^ are shown in FiglTJa). These numerical calculations differ by 
cutoffs Nb (the number of bands) and Nq (the number of order parameter plane- wave harmonics at different G) , and 
also illustrate the dependence on the orientation of q. By keeping the density at two fermions per site we invariably 
find that n(g) is linear for small q and reaches its minimum at a finite value oi q = \Q\ in the direction q ~ (1,1, l)/v3 
for the simple cubic cosine potential. It can be seen from ([3]) that by varying the interaction strength U all eigenvalues 
of the matrix 11 are shifted by the same amount. Pairing occurs when the smallest eigenvalue Ti{q) becomes negative, 
but this happens at a finite momentum Q. 

The momentum dependence of n(q) enters through both the denominator in ([3]) and the vertex functions. The 
former can lead to commensurate pairing instability if the interband pairing channels are strong enough, a particle- 
particle analogue of Halperin-Rice exciton condensation in indirect gap semiconductors^. The vertex functions, 
however, contain linear terms in q which respect lattice symmetries (such as qni, etc.), arising from phase-space 
restrictions for intraband and interband pairing. Their contributions are not canceled out in (|3]) thanks to the non- 
linear factor involving ^nk, and make Ii{q) ~ ^ql?! for small q leading to an incom,m,ensurate pairing instability. A 
negative value Aq < is required in order to secure a local minimum at finite Q. We never found a positive Aq in 
our calculations, which would imply a uniform pairing instability with an unphysical dispersion for Goldstone modes. 

Fluctuation corrections can be systematically included by promoting the theory to a simplectic Sp(2iV) symmetry 
group and performing a perturbative 1/iV expansion^. The corrections still involve the same vertex functions, so the 
linear-q features are unlikely to be accidentally canceled. Therefore, we can expect incommensurate finite-momentum 
pairing instability to all orders of \/N . 

However, an equilibrium superfluid phase with incommensurate pairing is not necessarily stable. The equilibrium 
order parameter must be a standing wave involving all symmetry-related wavevectors Q at which pairing occurs most 
readily. The obtained state is called a pair-density-wave (PDW) and breaks translational symmetry both with its 
global U(l) phase and density. An incommensurate density modulation is very frustrated in a lattice potential and 
thus not pinned to the lattice. Quantum fluctuations can remove this frustration by stabilizing some commensurate 
order, at least at Q = 0. The arguments above suggest that the linear-q dependence of n(q) can hardly be completely 
removed, but fluctuations can modify n(g) to push its quadratically-dispersing minimum to a finite commensurate 
wavevector. 



III. FLUCTUATION EFFECTS 



The emergence of a fluctuation-dominated regime can be clearly seen in renormalization group (RG) analysis^. The 
theory ([T]) contains a number of interacting fixed points associated with a diverging "vacuum" scattering length a at 



zero temperature. Any band-insulator can be viewed as a "vacuum" of particles or holes and taken to the unitarity 
limit by adjusting the chemical potential to its nearest band-edge and tuning U to a critical value U* appropriate 
for the given lattice potential. The universal theory near unitarity generally contains a single fermion species, either 
particles or holes, whose low-energy dynamics is entirely captured by an effective mass; fermionic excitations in other 
bands are gapped and thus irrelevant near the fixed point: 

S' = Jm fl^ {-^UJ + E{k)) 4 „ + uJvhVk,Vq flj,,+,jlj,,^,,p , 

where k = (w,fc), Vk — dLud'^k/(2Tr)'^^^, and E{k) ~ Eq + ^^. Since the only fermion field / in this theory 
lives in a vacuum state, the fermion propagator acquires no renormalization at all orders of perturbation theory (all 
Feynman loop diagrams are zero). Only the interaction vertex is renormalized by high energy fluctuations, from an 
exactly summable geometric series of ladder diagrams. This allows deriving the exact RG equations to all orders of 
perturbation theory^: 

^^2E, , ^ = (2 - d)C/ - ^C/^ (5) 

where Eg = Eq + U is the effective band-gap and 11 is a positive cutoff-dependent constant. The unitarity fixed point 
is E* — 0, U* — — ell^^ ind = 2-|-e>2 dimensions. Having the exact RG equations, we find the flow of couplings 
with the scale parameter: 




[i-i-ni7(o)]e'-nc/(o) 

In two dimensions any initially negative U flows toward — oo, while in three dimensions only sufficiently large 
attractive U undergoes a run-away flow. In both cases, however, the run-away flow escalates at a finite I where 
U{1) — > — oo has a vertical asymptote. This indicates a breakdown of the RG. The working assumption behind the 
RG was that high-energy fermions are not paired. When U — >■ — oo, even the high-energy fermions must form Cooper 
pairs, and this happens at a finite scale I. At very large attractive U one must abandon the theory ^ and continue 
the RG using a theory which contains bosonic fields at all energy scales. This clearly affects the universality class of 
the superconducting transition, as order parameter fluctuations become qualitatively important. 

The bosonic effective theory in question can be derived by rewriting the model ([1]) in the basis of Wannier states. 
Assuming that fermions live only at "high energies" , either due to their band-gap or pairing, we can integrate them 
out and arrive at a lattice boson model: 



ij ijkl 



(7) 



Here Kij(x) is a function analytic at x = 0, and the bosonic fields bi live on lattice sites and have multiple orbital 
states, both labeled by indices i,j. On the other side of the superfluid transition is a bosonic Mott insulator with 
gapped but long-lived bosonic excitations. 

This theory is valid up to the energy scales at which the boson excitations cannot decay into free fermion pairs. A 
bosonic insulator can be accommodated as long as the gap for collective excitations Af, lies below the threshold for two 
fermion creation 2Af. This threshold can be estimated in the mean- field theory from the frequency dependence of the 
bosonic Green's function G{q,il) = —Il~^{q,il.). The Figlljb) illustrates a Bose-insulator regime in two dimensions 
at T = in which the lowest excitations are infinite-lifetime gapped bosonic modes (darker shades correspond to 
larger positive 2Af ~ A^). The superconducting transition occurring from this bosonic insulator is fluctuation-driven. 
The transition is effectively in the BEG regime, even though the required attractive interaction U can be very small. 

The most familiar Mott insulators have an integer number of bosons localized at each site due to strong repulsive 
forces among them^. However, Mott insulators can be stabilized at other commensurate boson densities leading 
to translational symmetry breaking-i£. In the present model the density of bosons is an integer per site since they 
originate from pairing of fermions in a fully populated band. However, the model is inescapably a multi-orbital one, 
and there is no constraint on the density of bosons in any particular orbital. The multi-orbital character of the 
effective theory ([7]) is tied to the same properties of the vertex functions Q which led to the PDW pairing instability. 
Therefore, we conjecture that the microscopic models such as ([TJ and the corresponding effective theories ([7]) can 
contain stable bosonic Mott insulating phases, adjacent to the PDW superfluid, which break translational symmetry 
with density and orbital orders. 
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